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ABSTRACT 

Using the IRAM 30 m telescope, a mapping survey in optically thick and 
thin lines was performed towards 46 high mass star-forming regions. The sample 
includes UC Hll precursors and UC Hll regions. Seventeen sources are found 
to show "blue profiles", the expected signature of collapsing cores. The excess 
of sources with blue over red profiles ([A^biue ^ A^red]/Atotai) is 29% in the HCO"*" 
J=l-0 line, with a probability of 0.6% that this is caused by random fluctuations. 
UC Hll regions show a higher excess (58%) than UC Hll precursors (17%), indi- 
cating that material is still accreted after the onset of the UC Hll phase. Similar 
differences in the excess of blue profiles as a function of evolutionary state are 
not observed in low mass star-forming regions. Thus, if confirmed for high mass 
star-forming sites, this would point at a fundamental difference between low- and 
high-mass star formation. Possible explanations are inadequate thermalization, 
stronger influence of outflows in massive early cores, larger gas reserves around 
massive stellar objects or different trigger mechanisms between low- and high- 
mass star formation. 
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Introduction 



Inflow motion is a fundamental plienomenon during stellar formation. Although the 
search for inflow is usually more difficult than that for outflow, studies of inflow have made 
great progress since the 1990s. In low-mass star forming regions, inflow motions w ere detected 



at different evolutionary st ages, including Class -I, Class and Class I cores (IZhou et al. 



I993I : iMardones et~aDll997l : llee. Mvers. fc Tafalla 1^9991 : ICregersen et~aDl20od : Isvanslboosh . 
Recently, a number of inflow candidates were found in high mass star formation regions. 
Among a sample of 28 massive cores, 12 were found to show line profiles that peak at blue- 
shifted velocities (h e reafter "blue profiles"; see Sec t . 3.1) , the expected signature of infiow 
dWu fc Evan]l2003h . IPuUer. WiUiams. fc SridharanI J2005h (hereafter FWS05) detected such 
asymmetric profiles in 22 cores within a samp le of 77 high-mass proto-stellar objects (HM- 
POs). Most recently, IWyrowski et al.l (120061 ) detected 9 sources with a blue profile in a 
sample of 12 ultracompact (UC) Hll regions. 

Variation of infiow motion with time is critical for high mass star formation. It has 
been indicated that when a protostar reaches >10 Mq it can generate enough radiatio n 
pressure to halt spherical infall and inhibit its mass increasing (jWolfire fc Cassinellil 119871 ). 
Observationally, however, it is not yet clear how infiow is related to the evolution of massive 
(proto) stars. To study this problem, we have carried out a survey for a sample including 
both cores of UC Hll regions and precursors of UC Hll regions. 

While previous surveys using single point observations provided some statistical evidence 
for the occurrence of infall within massive cores, blue profiles can also be caused by rotation. 
Therefore maps of the molecular environment are indispensable. Mapping also allows us to 
locate the center of the infiow and to identify cores that are simultaneously showing evidence 
for in- and outfiow. 

Therefore, we conducted a mapping survey including 46 high mass star-forming regions 
which were selected applying three criteria: (1) The sources must have been mapped in 
the submillimeter or millimeter wavelengths with continuum or spectroscopy; (2) signal-to- 
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2OO2I ). With respect to their stellar content, we can divide the sample 
into two different groups of targets: (1) Thirty t hree sources lack 6 cm continuum emissio n 
and are precursors of UC Hll regions or HMPOs (IMolinari et al.ll2000l : iBeuther et al.ll2002l ). 
Among these, thirty are hosting a luminous IRAS source. The remaining three are associ- 
ated with IRAC (the InfraRed Array Camera on the Spitzer Space Telescope) point sources 
(W3-W and W3-SE) or are not hosting an IRAC source (18454-3). All 33 cores comprise 
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'group r. (2) Thirteen UC Hll regions are assigned to 'group IT. This letter presents a hst 
of the identified collapse candidates and provides the statistics of blue excesses. Detailed 
properties of individual cores will be analyzed in a future paper. 



2. Observations 

The observations were performed with the IRAM 30 m telescope at Pico Veleta, Spain, 
from July 28 to Aug. 1, 2005. Four receivers were used simultaneously, usually two at 
A~3 mm and two at A~1.3 mm (A/B configuration). For some sources, none of the four 3 and 
1.3 mm lines were optically thin. In these cases the tracer lines were changed employing two 
receivers at A~2mm and the other two to cover the upper part of the 1.3 mm window (C/D 
configuration). The lines and corresponding beam sizes, efficiencies, and channel widths are 
given in Table 1. The channel spacing and the bandwidth are 78.125 kHz and 105 MHz 
respectively. The weather was extremely good for summer conditions, allowing us to observe 
the HCO"^ J=3-2 transition at 268 GHz and leading to 3 and 1.3 mm (2 and 1.2 mm) system 
temperatures of order 150 and 400 K (200 and 550 K) on a scale. Pointing and calibration 
were checked by continuum measurements of the standard sources W3(0H), G34.24, and 
NGC 7027 and were found to be better than 4" and ±20%, respectively. All observations 
were carried out in a position switching mode. For each source we observed a nine point 
map in a cross pattern with a spacing of 15". If infiow signature was detected, the map 
was enlarged in most cases to cover the entire region showing this signature. The on-source 
integration time per position was 1 minute, yielding a Ta* lo" noise level of 0.07 K for the 
3 mm N^H + (1-0) line. For the data analysis, the GILDAS software package (CLASS/GREG) 



was used ( iGuilloteau fc Lucasll2000l ). 



3. Results and discussion 

3.1. Blue profile identification 

For self-absorbed optically thick lines, the classical signature of infiow is a double peaked 
profile with the blue-shifted peak being stronger, or a line asymmetry with the peak skewed 
to the blue side. While optically thin lines should show a single velocity component peaking 
at the line center. 

Among the 46 cores observed, five (054904-2658, G31.41+0.31, 18454-3, 18454-4, 19266-1-1745) 
will be ignored because they show either too complex spectral profiles, inhibiting a detailed 
analysis, or a lack of optically thin lines. Estimates of optical depths were obtained from 
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line ratios between different isotopomers of CO and CS and from tlie relative intensities 
of individual hyperfine components in the case of C^^O and N2H+. C^^O, C^''0, C^^S and 
N2H+ tend to be optically thin, while CS is optically thick. HCO+ opacities could not be 
estimated . However, the si r ailarit y of HCO+ and CS line shapes (see Sect. 3.2) as well as the 
results of iGregersen et al.l (I2OOOI ) and FWS05 clearly indicate that HCO+ is also optically 
thick. 

The 41 remaining sources were detected in at least one optically thick and one optically 
thin line. A blue profile caused by inflow motion with velocity v oc r~^/^ in a region with 
higher excitation temp erature (Tp..^) insid e requires Ta*(B)/Ta*(R) > 1. Here r is the radius 
of the collapsing core (IZhou et al.l Il993l ). Ta*(B) and Ta*(R) are the blue and red peak 
intensitie s of the optically thick line. We also define a dimensionless asymmetry parameter 
following iMardones et al.l (Il997l ). 5V = (Vthick-Vthin)/AV[;hin- VJ^hick is the peak velocity of the 
opaque line, V^hin and AT^hin denote the peak velocity and width of the optically thin line. 
Only for dV < —0.25 or > 0.25 the line profile is rated blue or red, respectively. 

Our sources (Table 2) discriminate among five main types of line shapes: (1) cores 
with lines showing a "blue profile" (in the following denoted with B); (2) cores with lines 
showing a "red profile" (R); (3) cores exhibiting blue and red profiles at different spatial 
positions (BRS); (4) cores where some lines show a blue profile, while others display a red 
profile (BRL); (5) cores without obvious asymmetric lines (S). Only cores showing at least 
one line of type B, but no lines of type R are identified as targets potentially undergoing 
infiow motion. 



3.2. Collapse candidates and their profile "excess" 

With the criteria outlined in Sect. 3.1, seventeen infiow candidates are identified (see 
Table 2). Ten belong to group I and seven are part of group II. To provide a typical example. 
Fig. 1 shows the infall signature of the group I core W3-SE. Fig. la displays the HCO^ (1^0) 
spectra, showing the angular size of the core. Fig. lb shows a number of profiles towards the 
central position. The HCO+ (1^0) and (3-2) lines as well as the CS (3-2) transition show the 
blue asymmetry. For the HCO"^ (1^0) line this is also demonstrated in the position-velocity 
(P-V) diagram of Fig. Ic. For comparison. Fig. Id shows a P-V diagram of the optically thin 
C^^O (1-0) emission. 



The quantity "excess" as defined by IMardones et al.l (119971 ) is E = {Nb — N^)/Nt, where 



A^B and Ar mark the numbers of sources with blue and red profiles. A^t is the total number 
of sources. For our survey the excess was calculated for the two HCO"*' transitions and the 
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CS(3-2) line. Fig. 2 shows the log[rA*(B)/TA*(R)] and 6V (see Sect. 3.1) distributions 
of the three individual lines. Statistical results are given in Table 3. The observed excess 
derived from the HCO"'" (1-0) and (3-2) lines is 0.29 and 0.11, respectively. Both are larger 
than those obtained by FWS05 for the same lines (0.15 and 0.04). For the CS transition we 
obtain 0.29. To evaluate the statistical significance of the determined values, we conducted 
the binomial test (see FWS05 and references therein). Probabilities that the excesses are a 
product of a random distribution are given in the last column of Table 3. These are 0.006 
and 0.01 for HCO+ (1-0) and CS (3-2) respectively. Apparently, both lines are sensitive 
tracers of potential inflow motion in massive cores. 

To evaluate differences between the two classes of cores (I and II; see Sect. 1) with 
respect to the excess, we used the HCO"^ (1^0) line, which was mapped in the largest number 
of sources. The results listed in the lower part of Table 3 include 16 sources with profiles of 
type B. The excesses observed for group I and II are 0.17 and 0.58, respectively. 

Twenty of our 46 sources overlap with those of FWS05. Among them are 19 group I 
sources (out of 33), but only one source is from group II (out of 13). Our study includes 
various CO and CS lines. We also made maps. Thus we can view the common objects from 
a different perspective and can check, how far the choice of different molecular transitions 
and the presence of maps is leading to contradictions with previously published results. 
Differences are indeed significant. For eight of the 19 overlapping type I cores we obtain 
different line asymmetry classification, emphasizing the need for detailed maps. Nevertheless, 
the overall difference in the HCO"*" (1-0) excess is negligible (0.17 versus 0.15). 

To summarize, both data sets indicate that the H CO"*" (1-0) exces s is low for UC Hll 



precursors. For UC Hll regions, our results and those of IWyrowski et al.l (120061 ) suggest that 
the excess is larger and more significant. From the binomial test for group I and II, the 
probability that the blue excesses (0.17 and 0.58 respectively) are arising by chance is 0.13 
and 0.008, respectively. 



3.3. A comparison with low mass star-forming surveys 



While low mass star-forming regions show infall from the Class -I to the Class I stages 
of evolution, high mass star-forming re gions also exhibit infall signatures from their earliest 



stages till a UC Hll region has formed (jWu et al.ll2005l : iBirkmann. Krause. fc Lemke 112006 



Oin et al.l 120061 ). In low mass cores the profile excess was found to b e 0.30, 0.31 and 0.31 for 



Class -I, and I core samples in the HCO"*" (3-2) line jEvansI (120031 ) and references therein) . 
There seem to be no significant differences among the cores in different evolutionary phases. 
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However, our samples show the excess of UC Hll regions far surpassing that of the UC Hll 
precursors. This may point to fundamental differences between low and high mass star- 
forming conditions. Possible causes to the higher blue excess in Group II sources may be: 
(1) The molecular gas surrounding UC Hll regions may be more adequately thermalized to 
show the blue excess, i.e. the excitation temperature of specific lines may increase more 
monotonically towards the center. Thus all lines may produce blue profiles indicating infall 
motion, while in younger cores still some lines may show red profiles. (2) The amount of 
dense cool gas is larger towards younger objects. Outflows of dense molecular gas may be 
more active around Group I objects, shaping more red profiles. (3) Low mass cores are 
relatively isolated and their gas suppl y is limited. Simulations showed that this may halt 



the increase of inflow (IVorobyovi l2005l ) . However, high mass stars form in giant molecular 
clouds and their inflow motions are not easily halted by the exhaustion of molecular gas 
before most of it is dispelled. (4) In low mass cores, star formation may be spontaneous. 
In high mass cores, collapse may be trigged by extrinsic disturbances and the collapse may 
take more time to develop. 

With respect to potential selection effects, we used the same criteria to identify the 
targets of the two separate groups of sources. Since this study is based on a limited number 
of sources, more data quantifying the blue excess as a function of evolutionary stage would 
be highly desirable. 



4. Outlook 

We have carried out a mapping survey towards 46 molecular cores associated with 
massive star formation. Seventeen collapse candidates were identified. Among them are 10 
UC Hll precursors and 7 UC Hll regions. Overall, statistical results indicate a predominance 
of blue over red profiles which is surprisingly similar to that obtained towards cores forming 
low mass stars. Among high mass star-forming sites, the probability to detect blue profiles 
seems to depend on evolutionary stage and increases from UC Hll precursors to UC Hll 
regions. Toward low mass star-forming sites, however, this effect is not observed, suggesting 
a more fundamental difference in the way stellar masses are assembled. Larger line surveys 
and more detailed maps in various molecular transitions are needed to improve statistical 
evidence in order to confirm or to reject this potentially important finding. 



We are grateful to the IRAM staff for their assistance and F. Wyrowski for useful 
discussions. This research is supported by the Grant 10128306 and 10733030 of NSFC. 
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Fig. 1.— W3-SE. a. HCO+ (1-0) grid. The X and Y axes denote R.A . and Dec. offsets in 
arcseconds relative to the NH3 peak position of iTieftrunk et al.l (119981 ). b. Spectra towards 
the central position. Molecular species and transition are given at the upper-right corner of 
each inset, c. Position-velocity diagram of the optically thick HCO^ (1^0) line along ADec. 
= 0. Contour levels are 0.5 K (3cr) and 1.0 to 5.8 K by 0.8 K. d. Position-velocity diagram 
of the optically thin C^^O (1-0) line along ADec. = 0. Contour levels are 0.3 K {3a) and 1.0 
to 2.5 K by 0.5 K. 
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Fig. 2. — The distribution of log [TX(B)/TX(R)] (left; ratio of blue versus red peak intensity) 
and SV (right; the line asymmetry parameter, see Sect. 3.1). The upper, middle and lower 
panels show results from HCO'^(l-O), HCO+(3-2) and CS(3-2), respectively. There is a 
total of 41 cores observed in the HCO+(1-0), 28 in the HCO+(3-2), and 28 in the CS(3-2) 
line. For the 6V histograms, BRS sources (for a definition, see Sect. 3.1) had to be excluded 
(6, 5, and 5 in the three transitions, respectively). There are additional 12, 8, and 10 cores, 
for which no blue versus red peak intensity ratio could be obtained. 
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Table 1. Observed lines^^ 



No. 


Line 


Frequency 


HPBW 










(GHz) 


(arcsec) 




(kms~^) 


1 


HCO+{l - 0) 


89.18852 


27.6 


0.77 


0.263 


2 


HCO+{3 - 2) 


267.55763 


9.2 


0.45 


0.112 


3 


CS{3 - 2) 


146.96905 


16.7 


0.69 


0.159 


4 


CS{5 - 4) 


244.93561 


10.0 


0.49 


0.122 


5 


N2H+{l-0) 


93.17378 


26.4 


0.77 


0.251 


6 


C^^Oil - 0) 


109.78218 


22.4 


0.75 


0.213 


7 


C^^O{2 - 1) 


219.56033 


11.2 


0.55 


0.137 


8 


C^^O{l - 0) 


112.35928 


21.9 


0.74 


0.209 


9 


^1^0(2- 1) 


224.71437 


10.9 


0.54 


0.133 


10 


C^^S{5 - 4) 


241.01618 


10.2 


0.50 


0.122 



a) HPBW: half power beamwitfi; rjiah- beam efficiency; 
AK-es^ cliannel width 
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Table 2. Profiles of the observed sources 



Source 








s'' 










Source 








6'' 




D= 






Name'* 


J2000 




J2000 


(kpc) Profile"* 


Ref. 


Name"* 


J2000 




J2000 


(kpc) Profile"* 


Ref. 


W3-W^ 


02 


25 


32 


4 


+62 


06 


01 


1.96 


B 


1 


18488+OOOOSE^ 


18 


51 


26 


6 


+00 


04 


07 


5.4 


BRL 


3,7 




02 


25 


39 


6 


+62 


05 


61 


2.3 


BRS 


1 


G34.26+0.16^^ 


18 


53 


18 


4 


+01 


14 


56 


3.7 


In, B 


6 


W3-SB^ 


02 


25 


54 


6 


+62 


04 


11 


2.3 


B 


1 


18621+0134^ 


18 


54 


40 


8 


+01 


38 


02 


5.0 


B 


3,7 


05358+3543^ 


05 


39 


10 


4 


+35 


45 


19 


1.8 


BRS 


3,7 


18530+0215' 


18 


55 


34 


2 


+02 


19 


08 


5.1 


S 


3,7 


05490+2658' 


05 


52 


12 


9 


+26 


59 


33 


2.1 




3,7 


S76E'-' 


18 


56 


11 





+07 


53 


28 


2.1 


S? 


4,9 


G10.47+0.03-" 


18 


08 


38 


2 


-19 


51 


50 


5.8 


B 


10 


18553+0414NE' 


18 


57 


53 


4 


+04 


18 


15 


0.6 


s 


3,7 


G12.42+0.50-" 


18 


10 


51 


8 


-17 


55 


56 


2.1 


B 


2 


19012+0536' 


19 


03 


45 


1 


+05 


40 


40 


4.6 


B 


3,7 


G12. 89+0.49" 


18 


11 


51 


3 


-17 


31 


29 


3.5 


BRL 


2,3,9 


19092 + 0841SW-' 


19 


11 


36 


7 


+08 


46 


20 


4.48 


BRL 


5 


013.87+0.28^^ 


18 


14 


35 


4 


-16 


45 


37 


4.4 


S7 


10 


G43.89-0.79^-' 


19 


14 


26 


2 


+09 


22 


34 


4.2 


B 


10 


18144-1723NW^ 


18 


17 


23 


8 


-17 


22 


09 


4.33 


R 


6 


19217+1651N^ 


19 


23 


68 


8 


+16 


57 


45 


10.5 


B 


3,7 


18182-1433^ 


18 


21 


07 


9 


-14 


31 


53 


4.5 


B 


3,7 


19266+1745^ 


19 


28 


64 





+17 


51 


56 


0.3 




3,7 


G19.61^^ 


18 


27 


37 


9 


-11 


56 


07 


4.0 


S7 


2 


19410+2336^ 


19 


43 


11 


4 


+23 


44 


06 


2.1 


B 


3,7 


18264-1152^ 


18 


29 


14 


3 


-11 


50 


26 


3.5 


R 


3,7 


S86SE^ 


19 


43 


49 


7 


+23 


28 


41 


1.9 


BRS 


4 


18306-0835^ 


18 


33 


21 


8 


-08 


33 


38 


4.9 


R 


3,7 


S87N^ 


19 


46 


20 


6 


+24 


36 


04 


2.3 


S? 


4 


G24.49-0.04" 


18 


36 


05 


3 


-07 


31 


23 


3.5 


B 


2,11 


20126+4104-' 


20 


14 


26 





+41 


13 


32 


1.7 


S7 


3,7 


18337-0743NE' 


18 


36 


40 


9 


-07 


39 


20 


4.0 


BRL 


3 


20216+4107' 


20 


23 


23 


8 


+41 


17 


40 


1.7 


S 


3.7 


18355-0650"'* 


18 


38 


14 


2 


-06 


47 


47 


4.2 


B 


8 


20319+3958' 


20 


33 


49 


3 


+40 


08 


45 


1.6 


S 


3,7 


18372-0541' 


18 


39 


56 





-05 


38 


49 


1.8 


B 


3,7 


22134+5834' 


22 


15 


09 


1 


+ 58 


49 


09 


2.6 


BRS 


3,7 


1S3S5-0512E' 


18 


41 


13 


3 


-05 


09 


06 


2.0 


R 


3,7 


23033+5951' 


23 


05 


25 


7 


+60 


08 


08 


3.5 


S? 


3,7 


G31.41+0.31-" 


18 


47 


34 


7 


-01 


12 


46 


7.9 




10 


NGC7538-11' 


23 


13 


44 


7 


+61 


26 


54 


2.8 


B 


2 


18454-3' 


18 


47 


55 


9 


-01 


53 


35 


5.6 




3,7 


NGC7538-N'^'^ 


23 


13 


45 


4 


+61 


28 


12 


2.8 


B 


2 


18454-4' 


18 


48 


01 


4 


-01 


52 


37 


5.6 




3,7 


23139+5939^ 


23 


16 


09 


3 


+69 


55 


23 


4.8 


BRS 


3,7 


18470-0044^ 


18 


49 


36 


7 


+00 


41 


06 


8.2 


R 


3,7 


23161+5912^ 


23 


17 


21 





+69 


28 


49 


5.7 


BRS? 


3,7 



^Indices "I" and "II", attached to the source names, denote group I and II, respectively (see Sect. 1). "*" denotes the optical thin 
line is quoted from Luo & Wu (2007, in preparation); b. The positions are taken from the references (last column), c. If there is 
no distance available for a core, the distance is calculated using the galactic rotation curve of Brand &: Blitz (1993). If there is a 

distance ambiguity, the nearer one is chosen, d. Type of detected line profile (Sect. 3.1 for difinitions) . References: 1. Ticftrunk et 
al. (1998); 2. Mueller ct al. (2002); 3. Bcuthcr ct al. (2002) and references therein; 4. Zinchenko ct al. (1997); 5. Molinari ct al. 
(2000); 6. Hunter et al. (1998); 7. Sridharan et al. (2002); 8. Wu et al. (2006) and the references therein; 9. Hughes &z. MacLeod 
(1994); 10. Hatchell et al. (2000); 11. Lockman 1989. 
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Table 3. Blue excess statistics 



Line/Source A^b Nt E p 

Results from selected lines 

HCO+ (1-0) 16 4 41 0.29 0.006 

HCO+ (3-2) 8 5 28 0.11 0.29 

CS (3-2) 9 1 28 0.29 0.01 

HCO"'"(1-0) line results w.r.t. source properties 

Group! 9 4 29 0.17 0.13 

Group II 7 12 0.58 0.008 

Total 16 4 41 0.29 0.006 



Nb- Number of "blue" profiles; iVR,: Number 
of "red" profiles; A^t: Total number of observed 
sources. E: Excess; p: statistical likelyhood that 
the result is caused by random fluctuations (see 
Sect. 3.2). Note that only 15 of the 17 inflow 
candidates show blue profiles in the HCO"*" (1-0) 
line. The other two are identified by CS (3-2) and 
HCO+ (3-2) spectrum, respectively. One of the 16 
HCO+(1-0) blue profils is a "BRL" source (see 
Sect. 3.1). 



